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Abstract 

This paper proves a version for stochastic differential equations of the Lie-Scheffers Theorem. 
This result characterizes the existence of nonlinear superposition rules for the general solution of 
those equations in terms of the involution properties of the distribution generated by the vector 
fields that define it. When stated in the particular case of standard deterministic systems, our 
main theorem improves various aspects of the classical Lie-Scheffers result. We show that the 
stochastic analog of the classical Lie-Scheffers systems can be reduced to the study of Lie group 
valued stochastic Lie-Scheffers systems; those systems, as well as those taking values in homogeneous 
spaces are studied in detail. The developments of the paper are illustrated with several examples. 



1 Introduction 

A differential equation is said to have a superposition rule (a more explicit definition is provided 
in the next section) whenever any of its solutions can be written as a given (in general nonlinear) 
function of the initial condition and of a fixed set of particular solutions. The first characterization 
of the existence of superposition rules was given by the Norwegian mathematician Sophus Lie in a 
remarkable piece of work [Lie93 where he established a link between the existence of superposition 
rules and what we nowadays call the Lie algebraic properties of the vector fields that define a time- 
dependent differential equation. This result is referred to as the Lie-Scheffers Theorem and systems 
that satisfy its hypotheses as Lie-Scheffers systems. 

Lie-Scheffers systems have been the subject of much attention due to their widespread occurrence 
in physics and mathematics. The reader is encouraged to check with [CGM00 [ ICGM07) . and refer- 
ences therein, for various presentations of the classical Lie-Scheffers Theorem, an excellent collection of 
examples of applications of this theorem, and for historical remarks. 

The main goal of this paper is the extension of the Lie-Scheffers Theorem to stochastic differential 
equations. This generalization is stated in Theorem l3.ll It is worth emphasizing that the main result of 
the paper, Theorem 13.11 cannot be seen just as a mere transcription of the deterministic Lie-Scheffers 
Theorem into the context of Stratonovich stochastic integration by using the so called Malliavin's 
Transfer Principle [Ma78 . This Principle states that whatever is true for standard differential equations 
also holds for Stratonovich stochastic differential equations; as we will see later on, there are purely 
stochastic conditions that appear in the statement of the theorem. 

Additionally, in proving Theorem 13. II we have carefully spelled out the regularity conditions needed 
for the result to be valid; those conditions are only vaguely evoked in the classical references or in the 
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cited papers that study the deterministic case. More importantly, a careful construction of the proof has 
lead us to realize that the hypotheses under which we can guarantee the existence of superposition rules 
can be weakened: the Lie algebra condition in the classical theorem can be replaced by an involutivity 
hypothesis that is, in general, less restrictive. 

The contents of the paper are structured as follows. Section [2] explains in detail the notion of 
superposition rule and includes a proposition that translates this concept into geometric terms. Section[3] 
contains the main theorem that we have already described. 

Section [4] is dedicated to the study of Lie-Scheffers systems on Lie groups and homogeneous spaces; 
this case is particularly relevant since, as we show in the first result of that section (Proposition 14. ljl . 
classical Lie-Scheffers systems (roughly speaking, those generated by vector fields that close a Lie alge- 
bra) can be locally reduced to this case via a theorem due to Palais. In that section we also show, as 
an example, how Levy stochastic processes can be seen as Lie group valued Lie-Scheffers systems. The 
section concludes with a brief presentation of the classical Wei-Norman method for solving Lie-Scheffers 
systems, adapted to the stochastic context. 

Section[5]contains a discussion on how the existence of a superposition rule for a stochastic differential 
equation makes available a remarkable feature that has deserved certain attention in the context of 
standard stochastic differential equations, namely, the fact that the stochastic flow can be written as a 
fixed deterministic function of the Brownian forcing of the equation in question. Indeed, a well know 
theorem by Ben Arous [B89] , that we state in the paper and whose proof is based on the use of stochastic 
Taylor expansions, shows that this property of the flow is available under exactly the same hypotheses as 
the classical Lie-Scheffers Theorem. Our main theorem allows, admittedly only to a certain extent, the 
generalization of this statement to any stochastic differential equation that satisfies its hypotheses; more 
specifically, any SDE generated by vector fields that span an involutive distribution has a superposition 
rule and hence its flow can be written as a fixed deterministic function of the initial conditions and of 
a set of solutions that contain the stochastic behavior of the resulting map. 

The paper concludes with a section that contains a number of examples that illustrate the develop- 
ments of the paper. 

2 Superposition rules for stochastic differential equations 

Let (fi, J 7 , P) be a probability space. We start by considering the stochastic differential equation 

ST = S (X, T) 8X, (2.1) 

where X : R + x ft — > R' is a given R'-valued semimartingale and S (x, z) : T x R l — > T p R™ is a 
Stratonovich operator from M. to R™. Sometimes we will choose a basis in T*M. 1 and will write down the 
Stratonovich operator S(x, z) in terms of its components (Si (x, z) , . . . , Si(x, z)) with respect to that 
basis. 

Definition 2.1 A superposition rule of the stochastic differential equation h2.1)) is a pair (<&, {Pi, . . . , 
r m }), where $ : R™( m+1 ) — ► R™ is a (not necessarily smooth) function and {T l : R + x £1 — > R" | i = 
I, . . . , m} is a set of particular solutions of \2.1\) such that any solution T of \2.1\) can be written, at 
least up to a sufficiently small stopping time r . as 

r = $(z 1 ,...,z™ ; r 1 ,...,r m ) =-. * (z-,r u . . . ,r m ) , 

where z — (z 1 , . . . , z n \ a set of n arbitrary constants associated with the initial condition of the solution 
T, that is, r(0,u;) = (z\ . . . ,z n ), for all u € f2. We extend to the stochastic context the terminology 
used for standard differential equations and we will call Lie-Scheffers systems the stochastic differential 
equations that admit a superposition rule. 



haza.ro and Ortega: Superposition rules and stochastic Lie-Scheffers systems 



3 



Remark 2.2 As we will see in examples later on in the paper, superposition rules exist only locally. 
That is why we can, without loss of generality, restrict our attention to stochastic differential equation on 
Euclidean spaces. Observe also that we are requiring that $ does not depend on time, the probability 
space, or the noise X. This prevents us from using certain regularization techniques at the time of 
testing the existence of superposition rules. For example, when dealing with a deterministic differential 
equation, the standard transformation of a time-dependent system 7 = / (t,-f) on R", / : R" ,+1 — > R n 
into the autonomous one 

7 = / (*) 7) an d t = 1 

on R™ +1 obtained by adding an extra trivial differential equation for the time, is not allowed; indeed, 
if we find a superposition rule for the transformed autonomous system, that rule does not yield a 
superposition rule for the original system that satisfies the requirements of our definition, precisely due 
to the explicit dependence on time that appears in the superposition function. 

In order to study the implications of the presence of a superposition rules we take a more geometric 
approach. Let \1/ be the function defined by 



(z,q a ,q 1 ,...,q m ) 1 — > q - $ (z; q u . 



(2.2) 



Notice that for any z G R™, the function * z := * (2, •) : R"( m + 1 ) R™ i s constant on a (to + l)-tuple 
(r,Fi . . . ,F m ) of solutions of the system (|2.ip . at least up to a given stopping time r, provided that 
r t= o = z G R n a.s.. From now on we assume that all the solutions T that we are dealing with are 
constant a.s. at t = 0. Additionally, if the function $ is smooth then the map ^ z : R n ( m+1 ) — > R" is a 
submersion for any fixed z G R™, because 



rank ( ? I = rank (J n ) = n (2-3) 



dql 



where /„ is the identity matrix of dimension n. Consequently, for any z 6 R™, the level set '3/ J 1 (0) C 
-§gn(m+i) j g a c i osec [ embedded submanifold of ]R"( m + 1 ) of dimension nm . That is, the function \1/ defines 
a family Q of regular nm-dimensional submanifolds Q z via the zero level sets ^J 1 (0) = {p G R n ( m+1 ) | 
"J (z,p) = 0} =: Q z of for any z G R n . The submanifolds Q z are globally diffeomorphic to R nm via 

the restriction 7r m | g ^ to Q z of the projection 7r m : R n ( m + 1 ) = R" x xl" — > R" m = R" x - m - x R" 
onto the last to 1™ factors. This is easy to see by verifying that the inverse 2 Z : R m ™ — > Q z of Tr m \g z is 
given by B. z (qi, ■ . ■ , q m ) = (^(2; qi, . ■ . , q m ), Qi, ■ ■ ■ , 9m), which is obviously a diffcomorphism. In order 
to study the significance of the family of submanifolds Q we start by introducing the following definition. 

Definition 2.3 Let Y : R™ -» R™ &e a vector field. The vector field 

(go, •••,5m) 1 — > {Y {qo) , • • • , Y (q m )) 
is called the diagonal extension ofY. 

It can be easily checked that the set of diagonal extensions of vector fields in X(R") are a subalgebra 
of X(R"(" i+1 )); more explicitly, for any Y U Y 2 , Y 3 € £(R") and A G R, 

[?i , % + AY3] = [Yi , Y2 + XY 3 ] . (2.4) 

The following proposition states that, roughly speaking, the family of submanifolds Q completely char- 
acterizes the superposition rule. 
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Proposition 2.4 Suppose that the stochastic differential equation \2.1\) admits a smooth superposition 
rule ($, {Ti, . . . , T m }). Suppose that (Ti, . . . , T m ) t= o — {p\, . . . ,p m ) G M. mn a.s.. Then, there exists a 
family Q of closed embedded nm- dimensional submanifolds o/K n ( m+1 - ) such that for any z G R n there 
exists Q z 6 Q such that (r z , T\, . . . , T m ) C Q z , with T z the solution of 12.1]) such that (T z ) t= o = z. 
Moreover, for any Q z G Q the map TT m \g z '■ Qz —> R nm is a diffeomorphism. 

Conversely, let Q be a family of (not necessarily embedded) submanifolds of K™( m + 1 ' diffcomor- 
phic to R' im via ir m and {Pi, . . . , T m } a set of distinct solutions of H2.1\) such that (F\, . . . , T m )t—Q = 
(pi, . . . ,p m ) G R" m a.s.. Then, if for any point z G K™ there is an element Q z that contains the point 
(z,pi, . . . ,p m ) and the diagonal extensions (S± (X, •),..., Si (X, •)) of the vector fields (S± (X, •),..., Si(X, •)) 
that define 12. 1\) are tangent to Q z when evaluated at (T z , Ti, . . . , T rn ), then \2. 1\) admits a (possibly non- 
smooth) superposition rule. 

Proof. In view of the remarks preceding Definition 12.31 we just need to prove that having a family Q 
that satisfies the hypotheses in the statement allows us to recover the superposition rule. 

Let -{Ti, . . . ,r m } be the set of fixed distinct solutions of ()2. 1|) . Denote pi = (Ti) t= o the (necessarily 
different) constant initial conditions of Tj, i = 1, . . . , m. Let z = (z 1 , . . . , z n ) G R™ be a point and let 
Q z be the submanifold in Q such that (z,pi, . . . ,p m ) G Qz] by hypothesis, this manifold is diffeomorphic 
to R ram via the map ip z = n m \g , where n m : R™( m+1 ) — > R™ m is the projection onto the last nm 
factors. In other words, the last nm coordinates of a point in K"0+i) serve as global coordinates of Q z . 
Introduce the projection 

: — ► R™ 
(q ,...,q m ) i — ► q a . 

We now define 

(T ) t M := tt°„ o ^ {(T 1 ) t (to),..., (r m ) t (w)) . (2.6) 

It is immediate to see that (ro) t=0 = z and that Tq is a semimartingale because, by construction, it is 
a composition of smooth functions with semimartingales. Let now T z be the unique solution of (|2.ip 
with a.s. initial condition z G 1". We will proceed by proving that Tq defined in (|2.6p equals T 2 and we 
will therefore have a superposition rule <E> given by the map $(z; Pi, . . . , T m ) := 7rjj„ o (pj 1 (ri, . . . , T m ). 
Notice that unless additional hypotheses are assumed on the family Q, there is no guarantee on the 
smoothness of $ on the z variable. 

In order to prove that Tq equals T z , denote by (q k ; k = 1, . . . , nj the coordinates on R" and by 
{ql; k = 1, . . . , n; a = 0, . . . , m) the coordinates on R n ( m + 1 ) . Let F£ : W lm -> E n and X% : R nm -> 
jjn(m+i) ^ ne ma p S defined as 

( d 

K • • • i *») = T (g 1 ,...,g m )( 7r R'' V* I'm) ( ^fc 

n entries 

X fe a (^(gi, . . . , q m )) = T iqi ,..., gm) (ip- 1 on m ) (J-^j = (F fc ° ( 9 i, . . . , q m ) , 0, »T. 1 ,(0,*7l ) l,...,0) ) m r.o, 0), 

where a = 1, . . . , m, k = 1, . . . , n. Observe that, by construction, the nm vector fields Xfi are linearly 
independent and span T q Q z at any q G Q z , since ipj 1 is a diffeomorphism form R 1 "™ to Q z . 
Now, we notice that for any j — 1, . . . , I, the vectors 

s i (x ; r" > r 1 ,... J r m ) = (5 i (Jf,r») >J s i (x,r 1 ) > ... > 5 i (x,r ro )) (2.7) 

are by hypothesis tangent to Additionally, due to (|2.6p and the Stratonovich differentiation rules 
we can write 

m n m n / 

= E E ^ a ^i. • • • > r ») ff ' = EEE^( r i- r ™) ^ r -) ( 2 - 8 ) 

a—lk—l a—lk—lj—1 
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Moreover, 

/ m n \ 

E E F k (Ti. ■ • • - r «0 5 " r -) . ^ (*> r i) . • • ■ - S i r m ) G R n(m+1) (2.9) 



\a=lfe=l / 

belongs also to TQ Z for any j = 1, . . . , Z, since (|2.9|) can be written as a linear combination of the nm 
linearly independent vector fields X£. Indeed, 

(m n \ m n 

E F k (Ti. • • • - r m) r„) , s,- (x, ro , . . . , Si (x, r m ) =EE 5 i ^- r «) x * ( r i> • • • > r ™) ■ 
a=lfe=l / a=lfc=l 

Subtracting (|2.9[) from (|2.7p , we see that for any j = 1, . . . , I, 

(m n \ 
(x,r*) - ^E F * (Tx, . . . ,r m )^ fe (x,r a ) ,o, . . . ,o e r&. 
o=lfc=l / 

Any of these vectors fields, if different from zero, is obviously linearly independent from all the X%, 
a = 1, . . . , m, k = 1, . . . , n. If that is the case we could therefore conclude that d\ra(Q z ) is strictly bigger 
than nm, which is obviously a contradiction. Therefore, Wj = necessarily, and hence 

m n 

Si (x, n = J2 E F k Fi, • • • ■ r m ) sj 5 (x, r„) , 

a=l fe=l 

which guarantees that To is a solution of (|2.1[) because by (|2.8p 

i 

<5r = ^ Sj (x, v z ) sx j = sr z . m 

3 = 1 

Remark 2.5 In the previous proposition we saw how the tangency of the diagonal extensions of the vec- 
tor fields that define the SDE to the submanifolds in Q is a sufficient condition to ensure the existence of 
a superposition rule. Is it necessary? Suppose that we have a smooth superposition rule ($, Ti, . . . , T rn ) 
and let be the associated map introduced in (|2.2|) . As we have that \l/ z (r z ,Ti, . . . , T m ) = 0, the 
Stratonovich differentiation rules yield 

i=l a =0 Ha j=l i=l a=0 Ha 

A sufficient condition for this identity to hold is that, for any j 6 {1, .., I}, 

n m o t 

^^^(f,^,...,^)^ (x, r„) = o (2.ii) 

i=l a=0 ?a 

or, equivalently, that the diagonal extensions Sj (X, T z , Ti, . . . , T m ) are tangent to the elements of the 
family of submanifolds Q given by the zero fibers of the maps \& z . Additionally, one can find situations 
in which ()2.10p implies (|2. 1 1[) : for instance if j — 1 and (like in the case of the Brownian motion) the 
quadratic variation [X, X] is a strictly increasing process, a straightforward application of the Doob- 
Meyer decomposition and the Ito isometry make in this case ()2.10j) and (|2.11|) equivalent. 
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Remark 2.6 If we add to the hypotheses of Proposition l2~4"l that for any z 6 K™ and for any (pi, . . . ,p m ) G 
R" m there exist a submanifold Q z in Q such that (z,pi, . . . ,p m ) G 5z (for instance when £7 is a foliation 
of R"( m + 1 ) whose leaves are diffeomorphic to 1™ via 7r m ) then the superposition function that we con- 
structed in the proof of that result has the following extremely convenient property: the superposition 
function is the same for any fundamental sets of solutions {Li, . . . , L m } that we may want to choose. 
In other words, once $ is know, we can take m arbitrary independent solutions of (|2.1[) to write down 
any solution. This situation frequently occurs in mechanics; see for instance, the study of the classical 
Riccati equation in ( |CMN98j ). 

3 The stochastic Lie-Scheffers Theorem 

The main goal of this section is proving a theorem that characterizes the existence of a superposition rule 
for a stochastic differential equation in terms of the integrability properties of the distribution spanned 
by the vector fields that define it. This can be translated into a Lie algebraic requirement, which allows 
us to recover the classical Lie-Scheffers Theorem in the stochastic context (Corollarv l3.4|) . 

In order to have at hand the necessary concepts to state the main theorem, we start by briefly recalling 
some standard results on generalized distributions due to Stefan |St74a[ ISt74bj and Sussman |Su73j . 
Let M be a smooth manifold, T> C X(M) be a family of smooth vector fields, and D the smooth 
generalized distribution spanned by T>. Let Gtj be the pseudogroup of transformations generated by 
the flows of the vector fields in V and constructed as follows: let k £ N* be a positive natural number, 
X an ordered family X = (X\, . . . , X)~) of k elements of £>, and T a fc -tuple T = (ti, . . . , i&) e R fc 
such that F\ denotes the (locally defined) flow of Xj, i G {1, . . . , k}, ti\ the elements Tt of Gtj are the 
locally defined diffeomorphisms of the form Tt = o o ■ ■ ■ o . Two points x and y in M are 
said to be Gx>-equivalent, if there exists a diffeomorphism Tt G Gtj such that Tt(x) = y. The relation 
Gt> -equivalent is an equivalence relation whose equivalence classes are called the Gx>-orbits, that are 
sometimes referred to as the accessible sets associated to the family T>. 

Given the family T> and the associated pseudogroup Gtj we can define another family £>' of vector 
fields as 

V := {TT T ■ X | X e V, T T G G v }, 

that clearly extends V, that is, T> C T>' . The distribution D' spanned by the elements of V is by 
construction Gtj -invariant. That is, for each Tt G Gtj and for each z G M in the domain of Tt, 

T z T t {D'{z)) = D'{T T {z)). (3.1) 

Moreover, since (£>')' = ^' by construction, the Stefan-Sussmann Theorem guarantees that it is com- 
pletely integrable in the sense that for every point z e M, there exists an integral manifold of D' 
everywhere of maximal dimension which contains z. The maximal integral manifolds of a completely 
integrable generalized distribution on M form a generalized foliation of M (see for instance [D85 ) . 
A leaf of a generalized foliation is regular if it has a neighborhood where the singular foliation induces a 
regular foliation by restriction. A point is regular if it belongs to a regular leaf. Regular points are open 
and dense in M ([E>85, Theoreme 2.2]). We will refer to D 1 (respectively V) as the Stefan-Sussmann 
extension of D (respectively T>) . The Stefan-Sussmann's Theorem also establishes an equivalence be- 
tween the Gp-invariance of D (£)' = D) and its complete integrability; additionally, if D is a completely 
integrable distribution, then its integral manifolds are the G£)-orbits. When the distribution D has 
constant dimension, the Stefan-Sussmann Theorem reduces to the celebrated and especially convenient 
Frobenius Theorem which states the D is integrable if and only if D is involutive. Recall that D is 
involutive if [X, Y] takes values in D whenever X and Y are vector fields with values in D. 
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In the sequel, we will use the following notation in order to be able to handle diagonal extensions 
of different dimensions. Given I £ N and X £ X(R"), we will denote by X 1 £ X(R'™) the diagonal 
extension of X to M. ln . For the sake of consistency with the previous section X means X m+1 . 

Theorem 3.1 (Lie-Scheffers' Theorem for SDE) Let 

ST = S(X,T)SX (3.2) 

be a stochastic differential equation on R", where X : R + x 57 — > R' is a given R' -valued semimartingale 

and S(x,z) : T X M. 1 > T p R" is a Stratonovich operator from R' to R n . Let V be an arbitrary open 

neighborhood of R" . Then, 

(i) // the X-dependent vector fields {Si (X, •),■■■, Si (X, •)} can be expressed on V as 

r 

Sj (X, z)=J2 b) (X) Y (z) e T Z W\ b) £ C°°(m. 1 ), zeV, (3.3) 

i=l 

and the distribution D spanned by the vector fields T> = {Yi, . . . , Y r } C X (V) is involutive, then 
\3.2\l admits a local superposition rule. 

(ii) Conversely, suppose that VS.Ufy admits a superposition rule (<&, {Fi, . . . ,r m }) and that the diago- 
nal extensions {Si (X, •),..., Si (X, •)} to R™( m+1 ) are tangent to the family Q of nm- dimensional 
submanifolds of R™( m + 1 ) associated to this superposition rule (see Proposition \2.4\ ). Let D(q) := 
span{Sj(X t , q) | j G {1,...,Z}, t 6 R+}, q £ M ri ( m + 1 ) ; D' the Stefan- Sussmann extension of 
D, and Go its associated generalized foliation. Let z £ R™, Pi — (Ti)t=o, and suppose that 
p = (z,pi, . . . ,p m ) £ R"( m + 1 ) belongs to a regular leaf (Go) z of Go. Then, there exists an open 
neighborhood V of z, a family of vector fields {Yi, . . . , Y r } C X (V), and a family of functions 

{ fi i>i=W c C °° ( R such that 

r 

S. ] (X,v) = Y J b){X)Y l (v), (3.4) 

for any v £ V . Moreover, the vector fields {Yi, . . . , Y r } form a real Lie algebra. 

Proof, (i) Given I £ N, we define V 1 := Vx. l ).xV ande?/ :— max qeV i |dim ^spanjY^g), . . . , Y^q)}^ |. 
Notice that for any I £ N one has di < di + i and di < r. Let m £ N be the smallest number for which 
d m = d m+ i and let go € V m+1 be such that 

dim (span^+^o), • • ■ , Y r m+1 (q )}) = d m+x . (3.5) 

The maximality of the dimension of span{Y 1 m+1 , . . . , Y" l+1 } at qo implies that there exists a neigh- 
borhood U of q in V m+1 for which dim(span{Y L m+1 (g), . . . , Y r m+1 {q)}) = d m +i, for all q € U. In- 
deed, the expression (|3.5[) is equivalent to saying that the r x n(m + 1) matrix M(q) with entries 
Mij(q) := (Y™ l+1 (q)y has rank d m when evaluated at qo which, in turn, amounts to the existence of a 
non- vanishing minor Md m+1 (qo) of M(qo) of order d rn+ i. Since the minor Md m+1 (q) depends smoothly 
on q and Md m+1 (qo) 7^ 0, there exists an open neighborhood U of go m V m+1 for which Md m+1 {q) ^ 0, 
for any q € U. This implies that dim(span{Y 1 m+1 (g), . . . , Y™ +1 (g)}) > d m +i, for all q £ U. How- 
ever, the maximality used in the definition of di+i implies that the previous inequality is necessarily an 
equality. 
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Consequently, we have found an open set U C V m+1 in which the distribution D spanned by 
the family {Y™ +1 , . . . , Y r m+1 } has constant rank. Moreover, l|2.4p and the hypothesis on {Yi, . . . , Y r } 
being in involution imply by the classical Frobenius Theorem that D is integrable. Let Go be the 
family of maximal integrable leaves of D that form a foliation of U m+1 . Now, shrinking U if necessary 
and using foliation coordinates for Go, wc extend the distribution D to another integrable distribution 
D D D of rank nm whose integrable leaves Q contain those of Qq, and for which the restrictions of 



pn(m+l) 



to the leaves in Q are diffeomorphisms onto their images. 



Let now {pi, ...,p m } be a set of m distinct points in V such that (pi,...,p m ) 6 n m (U) and 
{Fi, . . . ,T m } the solutions of of f|3.2[) such that (ri, . . . ,F m ) i=0 = (pi, ■ ■ ■ ,Pm) a.s.. Let F := (Ti, . . . ,T m ) 
and t the stopping time defined as r := inf{£ > | Tt ^ 7r m (t/)}. Since the vector fields 

r 

(x, r) = h ) ( x ) ^r +1 (r) 
j=i 

are tangent to the integral leaves of (?o an d hence to those of G, at least up to time r, Proposition 12.41 
guarantees the existence of a local superposition rule. 

(ii) We start the proof by providing a lemma that will be needed in our argument. 

Lemma 3.2 Let {Y±, . . . , Y r } C X(M") with r < mn and Zef {Yi, . . . , Y r } be the corresponding diagonal 
extensions to R"( m+1 ). Suppose that {T q TT m (Yi(q)) : . . . ,T q 7r m (Y r (q))} are linearly independent for any 
q in a neighborhood U C R™( m + 1 ). If the sum Y^i=i with b l £ C°° (U), i — 1, . . . ,r, is again a 
diagonal extension then the functions b % are necessarily the pull-back by u m of a family functions in 
C^i^m (U)). More specifically, if (q J a ;j = 1, . . . , n; a = 0, . . . , m) are coordinates for K"( m+1 ), then the 
functions {& 8 }j=i,.. r do not depend on (ql'tj = 1, . . . , ri). 



Proof. Using the coordinates (g : ' ;j = 1, 
!£{l,..,r},je{l, 



n}, such that the vector fields {Y±, 

d 



ri) for K n , there exists a family of functions A\ <E C°° (R n ) 
Y r } C X(R n ) can be written as 



Y(q) 



E 



A\{q) 



dqi 



which implies that the diagonal extensions have the expression 

m n 

Yi(qo,---,q m ) =^E^ 

Then, if we assume that 

]b % (q Q ,...,q m )Y(q Q ,...,q m ) = 



«=o.,=i 



d 
dq£' 



i=l 



i=i O =o j=x 



~A{{q a ) 



Oql 



is a diagonal extension on U, then there exist some functions {i? I }i=i,....r C C°° (R") such that 



E 



6 l (a ,...,g m )^(g Q ) = B»'( 9o ) 



r./ 



if 



a 



0, . . . , m, j 



1. 



,,n. 



That is, the r functions 6 J (go, • 



, q m ) solve the following subsystem of linear equations 
/ ^(qo, ■ ■ ■ ,q m ) 



( A(q ) \ 

A( gi ) 



( B(q ) \ 



V Aim) J 



V b r (q , 



\ B{q m ) j 



(3.6) 
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where A and B are the n(m + 1) x r and n{m + 1) x 1 matrices, respectively, defined as A(q a )ij = 
Aj(q a ) and B(q a )i = B l (q a ), a = 0, . . . , m. Now, the hypothesis on the linear independence of 

{Tn m (Yi), . . . , Tn m (Y r )} implies that the rank of the matrix (A(qi), . . . , A(q m )) is r < nm and hence (|3.6 
has a unique solution which coincides with the unique solution of the system 



/ A( qi ) 



V Air, 



V b r (q , 



( B(gi) 



V B(9m) 



(3.7) 



Since there is no dependence on the coordinates go in the augmented matrix associated to the sys- 
tem (|3.7|) . its solution (6 1 , . . . , b r ) does not therefore depend on go, as required. ▼ 



Suppose now that the stochastic differential equation (13. 2|) admits a superposition rule and that we 
are in the hypotheses of the theorem. We start by emphasizing that since the vector fields {Si (X, ■) , 
. . . ,Si (X, •)} are, by hypothesis, tangent to the elements of the family Q then their flows leave invariant 
those submanifolds and hence, the Stefan- Sussmann extension D' of D is also tangent to the elements 
of Q. This argument guarantees that, given the regular leaf (Go)z of Go, then there exists an element Q z 
in Q that contains it. 

Now since p — (z,pi, . . . ,p m ) G K™( m+1 ) belongs to a regular leaf (Go) z of Go, then there is an open 
neighborhood U of p where we can choose (taking regular foliation coordinates) a family of linearly 
independent vector fields {%., . . . ,%} C £(R ,l(m+1 )) that span the tangent spaces to the leaves of Go^U . 
The vector fields {Yi, . . . , Y r } can be chosen as the diagonal extensions of r vector fields {Yi, . . . , Y r } C 
X(R"), since the Stefan-Sussmann extension D' = spanjTJ^ • Si (X, ■) | i 6 {1, . . . , I}, Tt £ Gp} of D 
is made of diagonal extensions. Indeed, in order to see that D' is spanned by diagonal extensions, it 
suffices to notice that the flow F t of the diagonal extension Y £ X(R" ( ^ m+1 - ) ) of a vector field Y £ X(R") 
is F t (qo, ...,q m ) = (F t (go), • • • , F t {q m )), with F t the flow of Y; hence 

T q F t (Y(q)) - (T qo F t x . . . x T Qm F t ) (Y (q)) 

= (T qo F t (Y(q )), . . .,T qm F t (Y (g m )) = {TF^Y))(q) 

is again a diagonal extension. Given that by (|2.4p diagonal extensions form an algebra, the statement 
follows. 

Moreover, since the distribution D'\jj is regular and integrable then it is necessarily integrable in the 
sense of Frobenius, that is, there exist functions {c^}i.j y k=i,..,r C C' 00 (R™( m+1 )) such that 



i^l =E4-^- ( 3 - 8 ) 



fe=i 

Now, as [Yj, Yi] = [Yj, Yi], we conclude that X)I-=i c ji^k is a diagonal extension. Also, as the projection 
7r m is a local diffeomorphism when restricted to U f) Gz, the family of vectors {T7r m (Yi), . . . , TTT m (Y r )} 
is necessarily linearly independent. In these circumstances Lemma 13.21 implies that the coefficients 
i c ij}i ,j,k=i....r do not depend on the first n coordinates q^, j = 1, . . . , n. We now apply 7Tg„ (see (|2.5p ) 
on both sides of (|3.8p and we obtain 

r 

(v) = ■■•> Vm)Y k (v) (3.9) 

fc=i 

where « £ F := 7Tg„({7) and (gi, . . . , g m ) G R™ m is any arbitrary point such that (v, gi, . . . , g m ) G [7. 
Since the left hand side of (|3.9p does not depend on (gi, • • • , g m ) then the dependence of the coefficients 
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Cj^qi, . . . , q m ) on those coordinates is necessarily trivial which allows us to conclude that {Y\, . . . , Y r } 
close a Lie algebra. 

Finally, since the vector fields Sj (X, ■) are tangent to Goi j = then there is a family of 

X-dependent functions K (X, ■) G C°° (U) such that 

r 

S J (X } q) = Y / bj(X,q)Y l (q), 
i=l 

for any q G U. As Sj {X, ■) is also a diagonal extension, we can use again Lemma 13.21 in order to prove 
that the functions do not depend on go- Consequently, 

r 

% (X, q) = J2 h ) (91, • • • , 1m)) Y (p) . (3.10) 

i=i 

As we did in the previous paragraph, we apply 7r^„ on both sides of (|3.10[) 

r 

Sj (X, v)=J2 b) (X, . . . , q m )) Y (v) , 
i=i 

for any v S V. Again, we realize that since the left hand side of this equation is independent of 
(<7i, . . . , q m ), the dependence of the functions on the coordinates (q±, . . . , q m ) is necessarily trivial, 
which yields expression (|3.4p . ■ 

Remark 3.3 Theorem 13.11 is a generalization for stochastic differential equations of the classical Lie- 
Scheffers Theorem stated for time-dependent ordinary differential equations. That theorem claims that 
a differential equation y — Y (t, y) on K" given by a time-dependent vector field Y (t, •) e X (R n ), t € K, 
admits a superposition rule if and only if Y can be locally written in the form Y (t,y) = X)[=i /* (v) i 
where {f %=!, .. „ r C C° °(R) and {1^, . . . , F r } C X (M n ) form a (real) Lie subalgebra of (X(M), [•, •]) (see 
CGM07] and jCGMOOj ). In relation to the traditional presentation of the Lie-Scheffers Theorem, our 
Theorem O 

(i) weakens the hypotheses under which we can guarantee the existence of superposition rules. The 
involutivity of the vector fields {Yi, . . . , Y r } is, in general, less restrictive than requiring that they 
form a Lie algebra over the reals. We know a posteriori by the second part of Theorem 13 . 1 1 that . 
around regular points, if there exists a superpositon rule, the components {Si, . . . , Si} of the 
Stratonovich operator can also be expressed in terms of a family of vector fields that close a Lie 
algebra. 

(ii) carefully spells out the regularity conditions under which we have a converse; those conditions are 
only vaguely evoked in the already cited deterministic papers. 

(iii) It is worth noticing that, apart from the two points that we just explained, Theorem 13.11 cannot 
be seen as a mere transcription of the deterministic Lic-Schcffcrs Theorem into the context of 
Stratonovich stochastic integration by using the so called Malliavin's Transfer Principle |Ma78] due 
to the purely stochastic conditions that appear in the statement of the theorem. Those additional 
requirements have to do with the tangency of the diagonal extensions of the components of the 
Stratonovich operator to the family of submanifolds associated to the superposition rule (see also 
Remark [23]) . 
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In the next corollary, we show for the sake of completeness how the classical statement of the Lie- 
Scheffers Theorem (generalized to SDEs) can be easily obtained out of Theorem 13.11 

Corollary 3.4 Using the notation in Theorem \3. 11 suppose that the X -dependent family of vector fields 
{Si {X, •),..., Si (X, •)} that define the stochastic differential equation \2. 1\) can be expressed as 

r 

Sj (X, z) = J2 b ) {X)Y l {z)e T z R n , b) E C°°(M'), z E R" . 

i=l 

Let Lie{Yi, . . . , Y r } be the real Lie subalgebra of (X(R n ), [•, •]) generated by the family {Y±, . . . , Y r } C 
X(R"). IfIAe{Yi, . . . , Y r } is finite dimensional then \2.1\) has a superposition rule. 

Proof. Let D and D 2 be the generalized distributions associated to the families of vector fields 
V = {Yi, . . . , Y r } and V 2 = Lie{Yi, . . . , Y r }, respectively. Observe that if D{z) £ D 2 (z), z E M™, then 
since Lie{Yi, . . . , Y r } is finite dimensional, we can always complete the family {Yi, . . . , Y r } with a finite 
number of vectors {Zi, . . . , Z s } C T> such that D(z) = D 2 (z). We then write the X-dependent vector 
fields {S 1 (X, •),..., Si(X,-)} as 

r s 

Si {x, z) = J2 h ) ( x ) Y i W + E 4 W^)' z e Rn > 
»=i fe=i 

with aj = for any j = l,...,l and any k = l,...,s. Therefore, we may simply suppose that 
D(z) = span{Lie{Yi, . . . , Y r }(z)}, z E R" and since D 2 is trivially involutive, the corollary follows from 
Theorem O(i)- ■ 



4 Lie-Scheffers systems and stochastic differential equations on 
Lie groups and homogeneous spaces 

The Lie-Scheffers systems that are defined by a set of vector fields that generate a finite dimensional Lie 
algebra, that is, those that satisfy the hypothesis of Corollarv l3.4l or of Theorem l5.1l can be reformulated 
in the language of group actions. More specifically, as we see in the next proposition, such systems come 
down locally to studying the solutions of an equivalent Lie-Scheffers system on a Lie group. 

Proposition 4.1 Consider a stochastic differential equation that satisfies the hypotheses of Corol- 
lary \3.4\ Let z E M be a point such that there exists a neighborhood V of z in which the dimen- 
sion of Lie {Yi, . . . , Y r } is constant. Then, shrinking V if necessary, there exists a Lie group G such 
that dim (G) = dim (Lie{Yi, . . . , Yf}|y), a group action S : G x V — > V, and Lie algebra elements 
{£i > • • • j £r} C Q such that 

S(exp(^),z), zEV. (4.1) 

t=o 

Moreover, the solution starting at z E M of the restriction to V of the stochastic differential equation 
may be expressed as 

T z t = ~(g u z), (4.2) 
where gt : M + x ft — > G is the semimartingale solution of the stochastic differential equation on G 

r 

Sgt = J2^(9t)SXt (4.3) 

i=l 



with initial condition gt=o = e a.s. 
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Proof. Since the statement of the proposition is local we can always assume that the vector fields 
\Y\, .. . ,Y r } are complete by multiplying them by a compactly supported bump function and by re- 
stricting ourselves to an open neighborhood V consistent with that construction. In that situation and 
if dim(Lie{Yi, . . . ,F r }|y) < oo, Palais showed in |P57j (see Corollary in page 97 and Theorem III in 
page 95) that there exists a unique connected Lie group G contained in the group of diffcomorphisms 
of M and a left action 5 : G x M — > M such that (|4.1j) holds and T e E z : q — > Lie{Yi, ...,Y r }(z) is an 
isomorphism, for any z G V. 

Let now gt : K+ x O — > G be the solution semimartingale of the stochastic differential equation on G 



5g t = Yj?(g t )6Xl (4.4) 



i=i 



where £p 6 X (G) denotes the right invariant infinitesimal generator associated to £j G Q via the left 
translations of G on G. Given that any two infinitesimal generators £ G and £ M , £ G g, are related by 
the formula T g S z (£ G ) = £ M (3 (<?, z)), g G G, z G V, it is straightforward to verify that if g t is a solution 
of (|4.3p with initial condition gt=o = e a.s., then 



rf = s( 5t) z), 

is the solution of <5L t = 5^1=1 ^ (^t) such that r = z, 



Remark 4.2 Observe that (|4.2[) may be understood as a general reformulation of (|5.2p (see also |B89, 
Theoreme 19]). Processes of the type = S (g t , z) defined using a group action are sometimes called 
one point motions ( |L04j ) . 

The proposition that we just proved shows that for Lie-Scheffers systems defined by vector fields 
that generate a finite dimensional Lie algebra g, it is the associated Lie-Scheffers system on the Lie 
group G (|4.3|) that really matters. This is the subject of the rest of this section. 

Stochastic differential equations on Lie groups. Let now G be an arbitrary connected Lie group 
and g its Lie algebra. Let {£i, . . . , & } and {e 1 , . . . , e'} be dual bases of q and Q*, respectively. Left 
(respectively, right) translations on G will be denoted by L : G x G — > G (respectively, R : G x G — > G). 
With the same notation that we have used so far, let 

Sf> >5 ):T pfl ~ fl — T g G 

be a Stratonovich operator from g to G, where rj denotes the infinitesimal generator associated to the 
G-action on itself by left translations. Consider the stochastic differential equation associated to (|4.5[) . 

I 

Sgt = ^2^{gt)SXt, (4.6) 
i=i 

for some driving noise (semimartingale) X : M + x!l-t§. Using the equivariance of the vector fields 
£ G G X(G) with respect to right translations, that is, T^i? 9 (£ G (/i))) = £ G (R g (h)) for any g, h G G, and 
£ G q, it is immediate to check that if T e is the solution of (|4.6p with initial condition r|_ = e a.s., 
then the solution Tf starting at g G G is given by 



rf = L T .g = R g {T e t ) 



(4.7) 
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In other words, the stochastic differential equation (|4.6p has a superposition rule in the sense of Defini- 
tion \2.1\ and the superposition junction $ is given by 

$ : G x G — ► G 

(h,g) 1 — * L h g = R g h. 

It is also worth noticing that (|4.6[) is stochastically complete ( E82, Chapter VII §6]) since it is a 
left-invariant system. Therefore any solution of (|4.6p is defined for all (t, ui) G R+ x and, consequently, 
so is any one point motion and, in particular, any solution of any Lie-Scheffers system on a manifold M 
which can be globally considered as induced by a group action S : G x M — > M. 

Levy processes and Lie-Scheffers systems. This is an important class of Lie group valued stochastic 
processes and, as we will now see, a class of examples of Lie-Scheffers systems. Recall that a continuous 
process g : R+ x — > G is called a right Levy process if, for any = to < t\ < t% < ... < t n , the 
increments 

9t ,gt gt^,9ti9t^ ■ ■ -^tn-iSt" 1 ( 4 - 8 ) 

are independent and stationary. This means that the random variables in (|4.8[) are mutually independent 
and that their distributions only depend on the differences U — U—i, i G {1, . . . , n}. If g to ^ e a.s., we 
define gf — gtgZ. , which is a right Levy process starting at the identity. 

We are now going to see that continuous Levy processes and Lie-Scheffers systems are closely related. 
First of all, recall that any right Levy process on a locally compact topological group with a countable 
basis of open sets is a Markov process with a right invariant Feller transition semigroup {-Pt} tS R + given 
by P t f (g) := E [f (g%g)], g G G, where / : G — > R is any measurable function. Conversely, any right 
invariant continuous Markov process is a right Levy process ( [L041 Proposition 1.2]). Moreover, if g : 
R + x £1 — » G is a right Levy process, then there exists a Z-dimensional Brownian motion B : R + x £1 — > R' 
with respect to the natural filtration {^- t e } tg H + of the process g%, I = dim(g), with covariance matrix 
(ctij )% t j=i i and some constants {ci}j=i j such that 

/ (gt) = f (go) + E f t?[f] (g s ) SBl + ^ Ci f £ [/] ( 9s ) da, 

for any / G C 2 (G) and where, as before, {^i, . . . , £;} is a basis of g ( L04, Theorem 1.2]). This expression 
amounts to saying that the Levy process g : R + x il — > G satisfies the stochastic differential equation 

l i 

5g t = E (as) s s + E ep (,9s) 
i=i i=i 

and hence by Corollary 13.41 we can conclude that arty continuous right Levy process is a solution of a 
right invariant Lie-Scheffers system. Additionally, it can be shown in this context (see |L04[ Theorem 
1.2]) that one point motions obtained out of a G-action 5:GxM->M are Markov processes with 
Feller transition semigroup {Pt' I } t£R+ 

P t M f (z) = E [f(E (gt, *))] , zeM, feC (M) . 

Lie-Scheffers systems on homogeneous spaces. Let H C G be a closed subgroup of G and 
consider the homogeneous space G/H = {gH | g G G} with the unique smooth structure that makes 
the projection tth '■ G — > G/H into a submersion. The group G acts on G/H via the map A : GxG/H — > 
G/iJ on G/iJ defined by (h,gH) i— > (hg)H. It is immediate to check that the infinitesimal generators 
associated to the left G-actions on G and on G/if are 7r# -related, that is, 

T>h (e G (g))=? G/H (tth (g)) 
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for any jeG, any £ € q, and where £ G / ff (gH) — Jj| t _ A exp ( t ^) (ffif)- This straightforward observation 
has as an immediate consequence the next proposition: 

Proposition 4.3 Let X : R + x — > q be a g-valued semimartingale, G a Lie group, and H C G 
a closed subgroup. Let T be a solution of the Lie-Scheffers system defined by X and the Stratonovich 
operator |^.5| ) with initial condition r 4= o- Then, tth (T) is a solution of the Lie-Scheffers system on 
G/H 

i 

tr = J2t? /H Ft) 6X1 (4.9) 

3=1 

with initial condition tth (r t= o). 

Observe that since the Stratonovich operator (|4.5p is right invariant by the action of G, and therefore 
ii-invariant, and that since this action is free and proper, the previous proposition can be seen as a 
particular case of the Reduction Theorem in |LO07j . The next theorem is a transcription of the 
Reconstruction Theorem in LO07] into the present context and describes how to construct solutions 
in the opposite direction, that is, it tells us how to construct a solution T of the Lie-Scheffers system 
J4.6'[ j out of the solutions of two other dimensionally smaller Lie-Scheffers systems: first, a solution of 
the reduced system (j4.9[) and second, another solution of a new Lie-Scheffers system, now on H. 

Theorem 4.4 Let X : K + x SI — > g be a Q-valued semimartingale, G a Lie group, H C G a closed 
subgroup, and S the Stratonovich operator defined in J^.5[ ). Let R : H x G — > G be the (right) action of 
H on G by right translations and A an auxiliary principal connection on tth '■ G — > G/H . Then, any 
solution r of the system |^.6| ) can be written in the form 

r* = Rh t 9t = gtfh- 

In this statement, g : M + x SI — > G is a G-valued semimartingale horizontal with respect to A, i.e. 
J (A,Sgt) = G q, gt=a = ^t=a, and such that tth (<?t) is a solution of the reduced system |^.ff[ ). On 
the other hand, h : R + X Q — ► H is a H-valued semimartingale that satisfies the stochastic differential 
equation 

5h t = R(Y u h t )5Y t (4.10) 
with initial condition h t= o — e, and associated to the Stratonovich operator 

%,/i):T ? f, * T h H , , 

n T e R h {rf)=r 1 H {h), 1 ' > 

and the stochastic component Y : R + x!l-»| given by 

i 



2 — 1 

Proof. See |LO071 Theorem 3.2 and Proposition 3.4]. ■ 



4.1 The Wei-Norman method for solving stochastic Lie-Scheffers systems 

The method that we are going to develop in this subsection is a generalization to stochastic systems 
of the one proposed by Wei and Norman in [WN63, WN64 in order to solve by quadratures time 
evolution equations of the form — HtUt that appear in quantum mechanics, where both Ut and Ht 
are bounded linear operators on a suitable Hilbert space. This method has already been adapted by 
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Carinena and Ramos [CROlj to the study of deterministic Lie-Scheffers systems on Lie groups and it is 
their approach that we will follow. As we will see later on, the power of this method and the ease of its 
implementation depends strongly on the algebraic structure of the Lie algebra g of the group G where 
the solutions of the stochastic differential equation take values. 

Let r : M x £1 — > G be the solution of ()4.6j) such that T t= o = e 6 G a.s.; we write it down in terms 
of second kind canonical coordinates with respect to a basis • • ■ , 6} 01 the Lie algebra g. That is, 

r t = exp(4a)---exp(46), (4.12) 

where {d\, . . . ,d l t } is a family of real- valued semimartingales, d 1 : R+ x — > R, such that d\ =0 = 
a.s. for any i — 1, . . . , I. Notice that the expression (|4. 12[) is only valid up to the exit time of T from 
the neighborhood U e of e G G where the second kind canonical coordinates for G around the origin are 
valid. The key idea in this method is that if the functions d % were differentiable then 

(see [CR011 Eq. (33) and (34)]), where Ad 9 (?7) G Q is the adjoint representation of G on 9, g G G, -q G 9. 
In our setup we obviously cannot invoke the differentiability of the functions gP, however applying the 
Stratonovich differentiation rules to (|4. 12|) with d l our real- valued semimartingales, % = 1, . . . , I, we have 

ST t = T e R rt (^; =i 54 (n i<4 Ad exp(dkj) ) ft) • 

This expression implies that for any right invariant one-form /i G G fi(G), that is, ^ G {g) — T*R g -i(fj,) 
for any 5 G G and a fixed G g* , 

/ (^^L) = (M,E / (lL.An^v ^d|>. (4.13) 

i—l 

At the same time, it is clear that / (a* , ST) = (n, X) and hence (|4.13p implies that 

^=E/(n^ A d OXP( ^))«- 

Using the identity Ad oxp ( ?? ) = e ad ^'^ = J2n>o "h. ^(v) ■"■ ac K?7), for any 77 G Q, and writing A = 
Yli=i ^(.i, we get the relation 

i=l (=1 J 

The system of stochastic differential equations (|4.14l) can be solved for the semimartingales d\, i = 
1, . . . , m by quadratures if the Lie algebra g is solvable (sec [WN63 ( WN64 ) and, in particular, for 
nilpotent Lie algebras. The solvable case was extensively studied in TC80| where similar conclusions 
were presented using a different approach. 

As a simple example consider the affine group in one dimension _4i,that is, the group of affinc 
transformations of the real line. Any element of A\ can be expressed as a pair of real numbers (do, a\) 
with a\ ^ defining the affine transformation x <— * a%x + qq. The product * : A\ x A\ — » A\ in A\ is 



(a , ai) * (6 , h) = (ao + ai&o, ai^i) • 
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If {£o = (1, 0) , £i = (0, 1)} is a basis of the Lie algebra cii of Ai, it is immediate to check that 

= = (4.15) 

Furthermore, the infinitesimal generators associated to the left action of Ai on itself are 

A typical Lie system on Ai would be, for instance, the following Stratonovich differential equation on 
the upper half-plane H + = {(x, y) G R 2 \ y > 0}, 

ST x = dt + r x 8B u 5Y y = dt + T y 5B t 

obtained as a particular case of Q4.5JI when G = A\ and X = (t,B), where B : R + x — > R is a 
Brownian motion. More generally, let X : R + x f2 — > ai be an ai-valued semimartingale and write 
X = X°£ + X 1 ^, with X° and X 1 real semimartingales. Then, using (|4.15p . (|4.14|) reads in this 
particular case 



x% + = y + y (6 - h 1 = (y h° - y & + (y ^ 

Putting together the terms that go both with £i and £o respectively, we obtain 



6 



d\ = xl 4 = x?+ f 

Jo 



d°JXl 



and hence 



5 The flow of a stochastic Lie-SchefFers system 

Theorem 13.11 claims, roughly speaking, that the stochastic system (|2.1[) admits a superposition rule 
($, {Ti, . . . , T m }) if the components of the Stratonovich operator S(x,z) : T X M. 1 — > T p R n , x G M', 
p G R", that define it may be written as Sj(X,z) = £[ =1 6*. (X) Y t (z), where bj G C°°(]R') and 
{Yi,...,Y r } C X(R") span an involutive distribution. The converse of this statement is also true 
provided that, for a given initial condition z e R™, the point (z, (ri, . . . ,F m ) t=0 ) is a regular point of 
the foliation Qo generated by the diagonal extensions of {Si(X, ■),..., S m (X, •)}. Notice that this is a 
reasonable condition since the set of regular points of a generalized foliation is open and dense ( D85, 
Theoreme 2.2]). Moreover, when this happens, the vector fields {Y\, . . . , Y r } form a real Lie algebra. 

The condition on the vector fields {Yi, . . . , Y r } forming a real finite dimensional Lie algebra or, more 
generally dim (Lie{Yi, . . . , Y r }) < oo, are particularly appealing since these are algebraic requirements 
that we may expect to be easily verified for stochastic differential equations of a certain type. Moreover, 
these conditions have consequences that go beyond Corollarv l3.41 More especifically, we will show that if 
dim (Lie{Yi, . . . , Y r }) < oo, then the general solution of a stochastic differential equation can be written 
by composing a deterministic function with a suitable noise. In the following paragraphs we are going 
to give a precise meaning to this statement and to put it in the context of well known results available 
in the literature. 

Traditionally stochastic differential equations on a manifold M have been presented as 



r 

sr t = Yo(r t )dt + J2Y*(r t )SBi, 

i=l 



(5.1) 
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where {Yq, .. . ,Y r } C X (M) and B : R + x — > W is a r-dimensional Brownian motion denned on 
a standard filtered probability space (£l,J-t,P). For the sake of having a more compact notation, we 
writei?° :— t. The flow of such a stochastic differential equation may be locally written, that is, up 
to a given stopping time r, by means of a Taylor series expansion that comes out of Picard's iterative 
method for solving stochastic differential equations. In order to be more explicit we introduce some 
notation. Let J = {ji, . . . ,j n }, ji G {0, . . . , r}, 1 < i < n, be a multi-index of size n. \\J\\ will denote 
the degree of J that, by definition, is the size of J plus the number of zeros in the n-tuple (ji, . . . , j n ). 
For any J = . . . , j n }, we consider the iterated Stratonovich multiple integral 

B J t = /'••• /' SBil-.-SBt. 



0<ti<. ..<*„<* 



In addition, Yj will denote 



If y G X (M) is a vector field on the manifold M, we will use the following notation for its flow: 
exp (sY) (z) denotes the solution at time s of the ordinary differential equation 7 = Y(j) with initial 
condition 7(0) = z. Then, 

Theorem 5.1 (|B89, Theoreme 20]) With the notation introduced so far, if dim (Lie{Fo- • • • , Y r }) < 
00 and span{Lie{loi ■ ■ ■ >^r}} ^ as constant dimension on a neighborhood V of the point z £ M, then 
there exists a stopping time r such that the solution of \5.1\) with initial condition z can be expressed as 

rr = exp(f] w ( 5 - 2 ) 

\n=l ||J||=„ 

up to time r. In this expression, 0j :— ^2 aeS 2 f"-n ^tr(J}> S n denotes the permutation group of n 
elements, and e(a) is the cardinality of the set {j S {1, . . . , m — 1} | cr(j) > a(J + 1)}. 

If the finiteness condition on the dimensionality of the Lie algebra generated by the vector fields 
is not available but, nevertheless, {Yq, . . . ,Y r } are Lipschitz vector fields, then the solution of (|5.ip 
starting at z G M can always be approximated by a process like (|5.2[) : if denotes the finite sum 

E«=iE||j||=«/5j-B/ J then 

r^=e X p(cf) W + ^/^W 

where the error term Riy(t) is bounded in probability when t tends to QC93 : , Theorem 2.1]). The 
expression (|5.2p also holds if instead of the hypotheses of Theorem 15.11 we require M to be an analytic 
manifold and {Yo, . . . ,Y r } a family of real analytic vector fields (|B89, Theoreme 10]). An important 
consequence of Theorem 15.11 lies in the fact that the general solution of the stochastic differential 
equation (15. lj) may be written, at least locally and up to a suitable stopping time r, as the composition 
of a deterministic and smooth function, namely, the flow exponential, with the diffusion that defines 
the stochastic differential equation (sec H92J for a complementary reading). From this point of view, 
there is a strong resemblance between Theorem 15.11 and Theorem 13. II 

• First, by Corollary [33 all the systems that satisfy the hypotheses of Theorem 15.11 admit a super- 
position rule. 

• Second, the superposition rule allows us to write any solution as a composition of the deterministic 
function and the set of solutions {ri, . . . ,r m } that are responsible for the stochastic behavior of 
the resulting flow. 
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We conclude by quoting two references that study the nilpotent case (that is, the Lie algebra 
Liejloi ■ ■ ■ ,Y r } is nilpotent); this case has deserved special attention in the literature (see, for example, 
K80 ) because in that situation the Taylor series expansion of the flow in terms of iterated integrals 
in (|5.2j) becomes finite. We also recommend the excellent exposition in B04 for a complementary 
approach to the subject of Taylor series approximation of the general solution of (|5.1|k in this book it 
is shown that, for instance, the Carnot group of depth N = dim (Lie{Yo> ■ ■ ■ , Yr}) can be used in the 
nilpotent case to integrate the Lie algebra action of Lie{Y" , . . . , Y r } when one writes, as we did in the 
previous section, a Lie-Scheffers system as a stochastic differential equation on a Lie group that acts on 
the manifold in question. 

6 Examples. 

6.1 Inhomogeneous linear systems. 

Let Ak : R — » M„(R) a n x n time-dependent real matrix and Bk : R — > R™ a time-dependent vector 
for any k = 1, . . . , I. Let X : R + x il — > R' be a scmimartingale. An inhomogencuous linear system is a 
system of stochastic differential equations on R™ that may be written as 

i 

sr t = J2(MWt) - B k (t))sx* (6.i) 
fc=i 

Let (q 1 ,. . . , q n j be coordinates for R". It is an exercise to check that (|6.1|) can be equivalently written 

as 

In In 
k—1 i,j=l k—1 i,j=l 

where the vector fields Y^, Zj G X(R n ), i^j^k— 1, . . . , n, are given by 

Y l -o l — 7 ■- — 
' ~ q dqi ' ' ~ dqi ' 

Given that 

[Yi,Y l k ]=8 k j Y l i -5\Y^ [Y;,Z k ] = -SlZ J , and [Z t ,Z,]=0 

we see that the vectors {Yj,Z^ \ i,j,k = l,...,n} C X(R") span a Lie algebra isomorphic to the 
(n 2 + n)-dimensional Lie algebra of the group of affine transformations of R n . Therefore, the system 
f|6. 1(1 satisfies the hypotheses of Theorem 13.11 and hence it admits a superposition rule. In order to 
explicitly construct the superposition rule, let r e ^ be the solution of the homogeneous part of (|6.ip . 

i 

5T t = J2A k (t)(T t )6X k 

k=l 

with initial solution F^ = ej s R n a.s., where ej = (0, 0, 1, 0, . . . , 0) for any j = 1, . . . , n. Let T be 
a particular solution of (|6.ip with initial condition T t= o = G R n a.s.. Then, 

n 

r t = E^+F 4 

3=1 

is the general semimartingale solution of (|6.ip starting at z = (z 1 , . . . , z") 6 R n . 
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6.2 The stochastic exponential of a Lie group. 

Let G be a Lie group and q its Lie algebra. Let {£1, . . . , £/} a basis of q and X : R + x £1 — > g be a g- valued 
semimartingale. Observe that X can be written as X = Xa=i a t£i f° r a family of real semimartingales 
a 1 : R + xfl — > M, i = 1, . . . , I. Following HL86] and [EPOlj . we define the (left) stochastic exponential 
E(X) : R + x £1 — ► G of X as the unique solution of the Lie-Scheffers system on G given by 

i 

i=l 

with initial condition T t= o = e € G a.s.. Unlike the conventions used in Section 21 the vector fields 
(£i) G G X(G) here are not the right-invariant vector fields built from £j, i = 1, . . . , Z, but the left-invariant 
ones. That is, 

(&) G ( 5 ) = T e L g &), geG. 

Except for the fact that (£,i) G G 3-(G), i = 1, . . . , I, are now left-invariant, solving a Lie-Scheffers system 
on a Lie group such as those presented in Section [4] amounts to computing the stochastic exponential 
of a given g-valued semimartingale X. 

The stochastic exponential establishes a bijection between g-valued local martingales and martingales 
on G with respect to certain connections. Recall that, given an affine connection V : X(M) x X(M) — > 
X(M) on a manifold M, a M-valued semimartingale T : R + x £1 — ► M is said to be a V-martingale (or 
a martingale with respect to V) provided that 

f(T) - /(r t=0 ) - 1 /" Hess/(dF,dF) 

is a real local martingale for any / S C°°(M), where Hess / : X(M) x X(Af) — > C°°(M) is the bilinear 
form defined as 

He SS f(Y,Z) = Y[Z [f]]-VzY[f] 

for any Y, Z E X(M) (see |E891 Chapter IV]). When M = G is a Lie group, one can construct left 
invariant connections V by using bilinear skew-symmetric forms a : g x q — >K on the Lie algebra jj via 
the definition 

V c g?7 G := a(£,n), £,r) £ Q. 

The curves exp(t£) e G, where ^ S g and exp : g — > G is the Lie algebraic exponential, coincide with 
the geodesies c(£) with respect to these connections that start at e G G and that satisfy c(0) = ^. It can 
be shown ( |EP01|. Le mma 1.4]) that the connections built from a — and a(£, if) = \ [£,rj\ induce the 
same V-martingales on G. Moreover, with respect to these two connections, the set of V-martingales 
consists precisely of the processes of the form Tq£ (X) where X is a g-valued local martingale and 
Fo a G-valued J^-measurable random variable QEP01, Proposition 1.9]). This expresion provides the 
bijection between g-valued local martingales and V-martingales on G that we announced above. 

6.3 Geometric Brownian motion. 

Let (K+, •) be the Abelian Lie group of strictly positive real numbers endowed with the standard product. 
Its Lie algebra is simply R and, for any £ 6 R, the Lie algebra exponential coincides with the standard 
exponential, that is exp£ = e^; consequently, the infinitesimal generator (right or left-invariant) is 

£ K + ( q ) = for any q G R. 

Let G = R+ x .". x R + be the Lie group constructed as the direct product of n copies of (R+,-). 
Its product map ■ : G x G — > G is obviously (<ii, . . . , a n ) ■ (pi, . . . , b m ) — (ai&i, . . . , a n b n ), a,, bi G 
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R+ for any i = 1, . . . ,n, and its Lie algebra is g =TiR+ x .". x TiR + ~ R x .". x R = R". Let 
{£; = (0, 0, 1, 0, . . . , 0) | i = 1, . . . , n} be the canonical basis of g — R", \i = (fi 1 , . . . , [i n ), a = 
(a 1 , . . . , <r n ) 69a couple of elements of Q, B ; R + x Jl — > g a n-dimensional Brownian motion on some 
filtered probability space (f2, P, {J r t }t£K + ) , and consider the following Lie-Scheffers system on G 

5r t = U-V (L t )A + ^£f(r t )<5B t l , (6 .2) 
^ ' i=i 

where a 2 — ((cr 1 ) 2 , . . . , (<J n ) 2 ). Using coordinates (q 1 ,..., q n \ in G we can rewrite (|6.2p as 

tfg* - - ~(^) 2 ) qtdt + <J l ql5Bl, i = 1, . . . , n, 

which may be rewritten in terms of fto integrals as 

dq\ = fjqidt + (J l q\dBl, i=l,...,n. (6.3) 

The solutions of the n-dimensional system of stochastic differential equations (|6.3p are usually referred 
to as the geometric Brownian motion which is well-known for its use in the Black-Scholes theory of 
derivatives pricing as a model for the time evolution of the prices of n assets in a complete and arbitrage- 
free financial market. 

The well-known solution of the differential equation (|6.3p can be easily obtained by using the 
stochastic version of the Wei-Norman method that we introduced in Section 14.11 Indeed, let qt — 
exp(aj£i) • ■ • exp(a"£„) be the solution of (|6.3|) starting at e = (1, . . . , 1) e G as in the, where a 1 : 
R + x SI — > R are real semimartingales such that a\ =0 — a.s. for any i = 1, . . . , n. Since the Lie algebra 
g of G is Abelian, and (|6.2p is written in Lie-Scheffers form 

1 

»=i 

by taking the noise semimartingale X := (^M 1 — ^2 ) * + (jl ^t' • • • > (m™ — ^T^ - ) * + cr ™^™) > tnc 
equation (|4. 14[) in the Wei-Norman method reduces to 

n 

(M 1 - (a 1 ) 2 /2, . . . ,u. n - (a n f/2) t + (tP-Bl . . .,a n B?) = 

1=1 

which implies that a\ — (// — (a 1 ) 2 /2)t + <J % B\ for any i = 1, . . . , n. Now, since the exponential map is 
given by 

exp : g — > G = R'| 
i Vr,-, — (e^,...,e^) 

where e^ is the standard exponential function, we recover the well-known result that the general solution 
q t of (16. 3p starting at qo € R™ is 

ft = ( (7 ie^ 1 -( CTl ) 2 / 2 ) t+CTls s...,g"e^"~( CT ' , ) 2 / 2 ) t + CT " B " 
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6.4 Brownian motion on reductive homogeneous spaces and symmetric spaces. 

Let G a Lie group and H C G a closed subgroup. We say that the homogeneous space M = G/H is 
reductive if the Lie algebra g of G may be decomposed into as a direct sum g = t) © m where [) is 
the Lie algebra of H and m is a subspace invariant under the action of Ad/f . That is, Adh (m) C ra 
for any h G H and, consequently, [f),ra] C ra. Suppose now that the reductive homogeneous space M 
is Riemann manifold with Riemmanian metric 77 and that the transitive action of G leaves the metric 
77 invariant. We want to define Brownian motions on (M, 77) by reducing a suitable process defined on 
G. The notation and most of the results in this example, in addition to a comprehensive exposition 
on homogeneous spaces, can be found in [H78 and (KN69]. The reader is encouraged to check with 
ELL98] to learn more about the geometry of homogeneous spaces in the stochastic context. 
We start by recalling that a Af-valued process T is a Brownian motion whenever 

/(D-/(r )-i J A(/) (r s )d s 

is a real valued local semimartingale for any / G G°°(M), where A denotes the Laplacian. As is widely 
known, the Laplacian is defined as the trace of the Hessian associated to the Riemannian connection V 
of rj. That is, 

r 

A (/) (m) = (A* o C Yt ~ VyT.) (/)(m) 

i=l 

where {Yi, . . . , Y r } C X (M) is family or vector fields such that {Yi(to), . . . , Y r (m)} is an orthonormal 
basis of T m M, m G M. 

Let o G M denote the equivalent class of H in M. We have assumed that (M, 77) is a Riemann 
manifold wiht a (left) G-invarinat metric rj. Since 77 is G-invariant and <E> is transitive, the only thing 
that really matters as far as the characterization of 77 is concerned is the symmetric bilinear form 
770 : T M x T M — > T Q M . It can be easily proved that there is a natural one-to-one correspondence 
between the G-invariant Riemannian metrics 77 on M — G/H and the Ad#-invariant positive definite 
symmetric bilinear forms B on T Q M = g/f) ([KN69, Chapter X Proposition 3.1]). The correspondence 
is given by 

ri =B 2-^(6)), 

where £1, £2 € 0, n : G — * G/if is the canonical submersion, and £ M € X (M) denotes the infinitesimal 
generator associated to £ £ g. In addition, if M is reductive then the bilinear form B may be regarded 
as defined on m, B : m x m — > K, since T Q M is naturally isomorphic to m, which is an Ad#-invariant 
subspace of g. The Riemannian connection V of the metric 77 associated to such a bilinear form B is 
given by 

v^£ 2 M = i[£f,£ 2 M ]+(f/(6,6)) M , (6-4) 

f [KN69|. Chapter X Theorem 3.3]). In this expression £1 and £2 belong to m and U : m x m — > m is the 
bilinear mapping defined by 

2B (u (6,6) ,&)=b (a, [6,&] m ) + s (Ks,ei] m ,6) , 

where [•, -] m is such that [•, •] = [•, ■]. + [•, -] m with [•, ■]. e \) and [•, -] m G ra. A consequence of (|6.4I) is 
that the Laplacian A takes the expression A (/) (to) = X)i=i(^5 M ^ te>^i) M ) (/) ( m )> m e 

M = G/K, where {£f f , . . . is an orthonormal basis of T m M. 

The most important examples of reductive homogeneous spaces are symmetric spaces. In that case, 
G is the connected component of the isometric group I(M) C Diff(M) of the symmetric space (M, 77) 
containing e = Id. In order to identify the symmetric space (M, rf) with a reductive space, take o G M a 
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fixed point and let s a geodesic symmetry at o. Then the Lie group G acts on M transitively and, if H 
denotes the isotropy group of o, M is diffeomorphic to G/H ( |H781 Chapter IV Theorem 3.3]). Suppose 
that dim (G) < oo and let a : G — > G be the involutive automorphism of G defined by a (g) — s o & g o s 
for any g £ G, where $:GxM^G denotes as usual the left action of G on M. It is a matter of fact 
that T e o~ : g — > g induces an involutive automorphism of g. That is, T e o~ o T e a = Id but T e a ^ Id. Let f) 
and m be the the eigenspaces of g associated to the eigenvalues 1 and —1 of T e a respectively such that 
2 = E) © m. It can be checked that f) is a Lie subalgebra of g, 

ft,f)] C h, [f),m] C m, [m,m] C J), 

and Ad// (m) C m ([KN69, Chapter XI Proposition 2.1 and 2.2]). Morevoer, the symmetric space G/K 
has a unique affine connection V invariant under the action of G. This is actually the Riemannian 
connection ( [KN691 Chapter XI Theorem 3.3]) so that (|6.4p reads 

V ef C 2 M = 

for any pair of left-invariant vector fields (J^ 1 and ^ ■ 

Returning to the general case, let {£i, . . . , £ r } be a basis of m such that {T e ir (£i) . . . , T e ir (^ r )} is an 
orthonormal basis of T a (G/K) with respect to n Q and let {£f , . . . , C X (G) be now the corresponding 
family of right-invariant vector fields built from . . . , £ r }. Observe that {^f f (m), . . . ,^ 7 (m)} is an 
orthonormal basis of T m (G/K) due to the transitivity of the action and to the G-invariance of the 
metric rj. Consider now the Stratonovich stochastic differential equation 

r r 

S 9t = E £(9t)SBi + &) G (6.5) 

i=l i=l 

where (Bj, . . . , B[) is a E r -valued Brownian motion. The stochastic system (|6.5|) is by definition if- 
invariant with respect to the natural right action R : K x G — ^ G, (p) = gfc for any g £ G and 
k <E K. In addition, it is straightforward to check that the projection ir : G — > G/X send any right- 
invariant vector field £ £ X(G), ^ e g, to the infinitesimal generator £ M £ X (M) of the G-action 
$ : G x M M. Hence (|6.5j) projects to the stochastic system 

r r 

i=i i=i 

on M by Proposition 14.31 It is evident that the solutions of (|6.6p have as a generator the second order 
differential operator | ° ^-? M + ^ tei Ci) ) an( i they are therefore Brownian motions. 

Acknowledgments The authors are indebted to Jose F. Carinena and Eduardo Martinez for their 
valuable comments and suggestions. They acknowledge partial support from the French Agence National 
de la Recherche, contract number JC05-41465. J. -A. L.-C. acknowledges support from the Spanish 
Ministerio de Education y Ciencia grant number BES-2004-4914. He also acknowledges partial support 
from MEC grant BFM2006-10531 and Gobierno de Aragon grant DGA-grupos consolidados 225-206. 
J. -P. O. has been partially supported by a "Bonus Qualite Recherche" contract from the Universite de 
Franche-Comte. 

References 

[B04] Baudoin, F. [2004] An introduction to the geometry of stochastic flows. Imperial College Press. 



Lazaro and Ortega: Superposition rules and stochastic Lie-Scheffers systems 



23 



[B89] Ben Arous, G. [1989] Flots et series de Taylor stochastiques. Probab. Th. Rel. Fields 81, 
29-77. 

[CGMOO] Carincna, J. F., Grabowski, J., and Marmo, G. [2000] Lie-Scheffers Systems: A Geometric 
Approach. Napoli Series on Physics and Astrophysics. Vol. 3. 

[CGM07] Carincna, J. F., Grabowski, J., and Marmo, G. [2007] Superposition rules, Lie theorem, and 
partial differential equations. Reports on Mathematical Physics, 60 (2), 237-258. 

[CMN98] Carinena, J. F., Marmo, G., and Nasarre, J. [1998] Int. J. Mod. Phys. A 13, 3601-27. 

[CR01] Carinena, J. F. and Ramos, A. [2001] A new geometric approach to Lie systems and physical 
applications. 

[C93] Castell, F. [1993] Asymptotic expansion of stochastic flows. Probability Theory and Related 
Fields 96, 225-239. 

[D85] Dazord, P. [1985] Feuilletages a singularites. Nederl. Akad. Wetensch. Indag. Math., 47, 21-39. 

[E82] Elworthy, K. D. [1982] Stochastic differential equations on manifolds. London mathematical 
society. Lecture notes series, 70. Cambridge Univ. Press. 

[ELL98] Elworthy, K. D., Le Jan, Y., and Li, X.-M. [1999] On the geometry of diffusion operators and 
stochastic flows. Lecture Notes in Mathematics 1720. Springer. 

[E89] Emery, M. [1989] Stochastic Calculus in Manifolds. Springer- Verlag. 

[EP01] Estrade, A., Pontier M. [2001] Backward stochastic differential equations in a Lie group. 
Seminaire de probabilites (Strasbourg). Tome 35, 241-259. 

[HL86] Hakim-Dowek, M. and Lepingle D. [1986] L'exponcntielle stochastique des groupes de Lie. 

Seminaire de probabilites (Strasbourg). Tome 20, 352-374. Lecture Notes in Math., 1204, 
Springer Verlag. 

[H78] Helgason, S. [1978] Differential geometry, Lie groups, and symmetric spaces. Pure and applied 
mathematics. Academic Press. 

[H92] Hu, Y.-Z. [1992] Serie de Taylor stochastique et formule de Campbell-Hausdorff, d'apres Ben 
Arous. Seminaire de probabilies (Strasbourg). Tome 26, 579-586. 

[KN69] Kobayashi, S. and Nomizu, K. [1969] Foundations of differential geometry (vol. II). Tracts in 
mathematics 15. Interscience Publishers. 

[K80] Kunita, H. [1980] On the representation of solutions of stochastic differential equations. 
Seminaire de probabilites (Strasbourg). Tome 14, 282-304. 

[LO07] Lazaro-Carm, J. -A., Ortega, J. -P. [2007] Reduction, reconstruction, and skew-product decom- 
position of symmetric stochastic differential equations http://arxiv.org/abs/0705.3156. 

[L04] Liao, M. [2004] Levy processes in Lie groups. Cambridge tracts in mathematics, 162. Cam- 
bridge university press. 

[Lie93] Lie, S. [1893] Vorlesungen iiber Continuierliche Gruppen mit Gcomctrischcn und Andcren 
Andwcndungen. Edited and revised by G. Scheffers, Teubner, Leipzig. 



Lazaro and Ortega: Superposition rules and stochastic Lie-Scheffers systems 



24 



[Ma78] Malliavin, P. [1978] Geometric diffcrentielle stochastique. Seminaire de Mathematiques 
Superieures, 64. Presses de PUnivcrsitc dc Montreal. 

[N66] Nagano, T. [1966] Linear differential systems with singularities and an application to transitive 
Lie algebras. J. Math. Soc. Japan. Vol. 18 (4), 398-404. 

[P57] Palais, R. [1957] A global formulation of the Lie theory on transformation groups. Mem. Am. 
Math. Soc. 22, 95-97. 

[St74a] Stefan, P. [1974a], Accessibility and foliations with singularities. Bull. Amer. Math. Soc, 80, 
1142-1145. 

[St 74b] Stefan, P. [1974b], Accessible sets, orbits and foliations with singularities. Proc. Lond. Math. 
Soc, 29, 699-713. 

[Su73] Sussman, H. [1973], Orbits of families of vector fields and integrability of distributions. Trans. 
Amer. Math. Soc, 180, 171-188. 

[WN63] Wei, J. and Norman, E. [1963] J. Math. Phys. 4, 575-81. 

[WN64] Wei, J. and Norman, E. [1964] Proc Amer. Math. Soc. 15, 327-34 



